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FIELD OF THE INVENTION 

This invention relates to imaging and in particular to 
methods and systems for image enhancement, 

BACKGROUND OF THE INVENTION 

Imaging involves transfer from the object domain into the 
image domain, but owing to limiting factors such as the 
finite size of energy source, detector size, sampling 
frequency, display density, software filter function, and 
possibly partial-volume effects experienced with some 
imagers, an infinitely fine delta function in the object 
domain cannot be faithfully reproduced in the image 
domain. Instead, a smeared-out image, or point-spread 
function (PSF) , is observed. Similarly, an infinitely 
sharp edge-response function (ERF) in the object domain 
becomes a smeared-out ERF in the image domain. The 
smearing effect becomes more intense as the adjacent ERFs 
of discontinuities or contrast profiles get closer to 
each other . 

It is an object of the present invention to provide a 
method and system by which the above problem can be at 
least partly overcome. 

SUMMARY OF THE INVENTION 

According to one aspect of the present invention there is 
provided a method wherein a de-convolution process is 
applied to the image-domain results of an o"bject-scan to 
derive therefrom the respective point- or line-spread 
function effective in the object- to image-domain 
transfer of one or more object-discontinuities, and to 
derive from said function the location in the image 
domain of the respective discontinuity. . 



According to another aspect of the invention there is 
provided an imaging system comprising means for 
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performing a de-convolution process on the image-domain 
results of an object-scan to derive therefrom the 
respective point- or line-spread function effective in 
the object- to image-domain transfer of one or more 
5 object-discontinuities, and means to derive from said 
function the location in the image domain of the 
respective discontinuity . 

The method and system of the invention enable the 
10 location of the respective discontinuity in the image 

domain, to be established with a high degree of accuracy. 
This is critical to image definition free of any 
substantial smearing, and to this end the location of the 
respective discontinuity may be derived to sub-pixel 
15 accuracy simply from the mid-point of the full-width 
half -maximum of said function. 

The said function may be correlated with the image-domain 
results of said transfer for enhancement of spatial 

20 resolution of the imaging of the one or more 

discontinuities. This enhancement may involve transfer 
of sub-pixels within the image-domain results of the 
respective one or more discontinuities, the sub-pixels 
being transferred within their respective image-domain 

25 results from one side to the other of said location. 

The* de-convolution process may be carried out using 
least-squares running filtering. 
BRIEF DESCRIPTION OF THE DRAWINGS 
30 An imaging method and system according to the present 

invention will now be described, by way of example, with 
reference to the accompanying drawings, in which: 

Figure 1 illustrates schematically the method and system 
35 of the invention; 
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Figure 2 illustrates features of processing performed in 
the method and system of Figure 1; 

Figure 3 shows results achieved from use of the method 
and system of Figure 1; 

Figure 4 shows to an enlarged scale a section of the 
contour of an image profile depicted in Figure 2; 

Figures 5 and 6 are a plan view and sectional end- 
elevation of a couch-top used in the method and system of 
Figure 1; and 

Figure 7 provides illustrates of a convolution operation, 
as a basis for a mathematical model of de-convolution 
processing in accordance with the method and system of 
Figure 1. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The method and system to be described with reference to 
Figure 1 utilise MR scanning for medical diagnostic and 
treatment-planning purposes . In principle and in the 
general techniques described, the method and system of 
the invention can be used in other applications of MR 
scanning and also in circumstances where other scanning 
techniques are utilised. Furthermore, although both 
structure and function are represented by discrete 
x boxes' 1 to 19 in Figure 1, the method and system are to 
a substantial extent manifest in programmed digital data- 
processing operations. 

Referring to Figure 1, data derived in accordance with 
conventional operation of an MR scanner 1 is processed 
for imaging purposes within a processor 2 . The output of 
the processor 2 is used to provide a display 3, and from 
this is subject to post-imaging processing 4. The post- 
imaging processing 4 includes the facility for selecting 
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a region of the display 3 for more-detailed and closer 
inspection. 

To the extent the imaging method and system of Figure 1 
have so far been described, they are conventional, and it 
is in further processing 5 of the image data of the 
selected region of interest obtained by the post-imaging 
processing 4 that a step forward from what is already 
known is achieved. More particularly, the further 
processing 5 is operative to define more clearly the true 
edges or boundaries of image contour (s) in the selected 
region of interest, and to enhance the accuracy of the 
imaging of those contours. 

The definition and accuracy of transfer of features from 
the object domain scanned by the scanner 1, to the image 
domain manifest in the post-imaging processing 4, is 
limited by many factors. The limitations arise from 
within the scanner 1 itself (in particular the finite 
size of the energy source) , within the processing 
performed by the processor 2, and within the display 3; 
limitations arise inherently from, for example, the data 
sampling frequency and display density used, and also 
from the filter-function of the software involved. More 
particularly, and referring to Figure 2, an infinitely 
fine delta function 20 in the object domain is not 
faithfully reproduced in the image domain. Instead, the 
transfer as represented by the arrow 21 results in a 
point-spread function (PSF) or smeared-out image 22 in 
the image domain. Similarly, an infinitely sharp edge- 
response function (ERF) or step 2 3 in the object domain 
becomes through the transfer represented by arrow 24, a 
smeared-out transition 25 of spread represented by 
dimension arrows 26, in the image domain. When two image 
ERFs are close to one another, the smeared-out effects 
run into each other. The consequent deterioration of the 
spatial resolution is often monitored .by the percentage 
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modulation transfer which is given by the ratio, 
expressed as a percentage, of the amplitude of the 
modulation in the image domain to that in the object 
domain. 

The smearing effect becomes more intense as adjacent ERFs 
of discontinuities or contrast profiles get closer to 
each other (or as the spatial frequency of the modulation 
becomes higher) ; this also causes loss of profile height. 
The inherent loss of the spatial resolution (that is, the 
part that is indicated by the smeared-out effect on the 
corner edge of the ERF) cannot be restored or partially 
restored even by re-scanning the image with an ultra high 
resolution digital scanner system. 

The further processing 5 is operative in accordance with 
the invention to provide accurate image-edge definition 
and location, and to improve spatial resolution in the 
imaging. More especially, in the context of Figure 2, 
the edge position corresponding to the discontinuity or 
step 2 3 of the object ERF is pin-pointed in the image 
domain from the mid-point of the full-width half -maximum 
(FWHM) of the image PSF; the pin-pointing is to sub-pixel 
accuracy for the image ERF. Low-contrast and ^area' 
filtering are used to remove x spurious' edges, and sub- 
pixel sampling to detect detail to the resolution of the 
single-pixel modulation. The discontinuity or step 23 of 
the ERF is then restored within the image domain by 
removing the sub-pixel values from outside the optimum 
edge position to compensate for those within. It is to 
be noted that the sub-pixels then become pixels in 
display, and that the enhancement is equivalent to the 
performance of an extra high resolution image transfer 
system. 

As represented in Figure 2 by the arrows 27 and 28, the 
image ERF 25 of an infinitely-sharp step 23 can be 
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produced by convolution of the image PSF 22 with the 
object ERF 23. In accordance with the present invention, 
de-convolution of the image ERF 25 using sub-pixel 
sampling, represented by the arrow 29, reproduces the 
image PSF 22 in a de-convolution space as image PSF 22'. 
The image ERF 2 5 is superimposed on the PSF 22' within 
this space as image ERF 25', and the optimum edge- 
position 3 0 is derived from the mid-point of the FWHM of 
the image PSF 22', and is pin-pointed to sub-pixel 
accuracy . 

For one-dimensional cases, the operation in accordance 
with the invention is relatively simple, as only either 
the x- or the y-profile, that is to say a line spread 
function LSF is involved. But for two-dimensional 
operations, both the x- and y-profiles, and if necessary, 
the xy-diagonal profiles to eliminate any possible 
streakings in the image, may be used; in this case, a 
proper weighting scheme will be required to re-construct 
the image. 

Once the original sharp-edge feature represented by the 
object ERF 23 is pin-pointed at the position 3 0 within 
further processing 5, that feature may be restored by 
additional re-processing 6 (Figure 1) . In re-processing 
6, the sub-pixel values occurring x outside' the optimum 
edge-position 3 0 are transferred to compensate those 
x within'. This is illustrated in Figure 2 by arrow 31 
transferring sub-pixel blocks 32 from after point 30 in 
the image ERF 25', to before it. The re-construction of 
image ERF 25' into image ERF 3 3 conforming closely in 
configuration to object ERF 2 3 is represented by arrow 
34. Image ERF 3 3 is displayed in enlarged form in 
display 7 (Figure l) . 

These techniques enable substantial "recovery of the loss 
of spatial resolution in the imaging, without the 
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trade-off loss of other properties such as image noise. 
Furthermore, the enhancement of spatial resolution in 
display 7 reproduces the region of interest selected from 
display 3, without blurring' (or step) effects at the 
profile edge. 



Figure 3 is illustrative of some of the low-contrast 
results provided in practice from MR scanning of a pig's 
brain in fluid. Curve A is the image ERF produced, 
whereas curve B is the line spread function (LSF) 
resulting from de-convolution of curve A carried out in 
processing 5. The optimum edge-position is established 
from the mid-point C of the FWHM of curve B, and the 
additional re-processing 6 is operative by means of sub- 
pixel transfer, to re-construct curve A to conform 
substantially to the edge-feature from which it 
originated in the display 7. 

It is to be noted that whereas curve A is stepped, curve 
B is nonetheless smooth and that mid-point C is located 
to sub-pixel accuracy. Furthermore curve B indicates a 
sensitivity of more than 8:1 between the profile-height 
and background noise. 

Referring again to Figure 2., the complete profile 3 5 of 
an image within the selected area of interest of display 
3, is built up as indicated by arrow 36, from the edge- 
position data derived within processing 5. This data 
identifies the location of the point 30, together with 
the locations of all corresponding points derived from 
sampling the multiple x- or y-profiles of the selected 
area of interest. The build up and display of these . 
points from the data takes place in display 8 so that a 
substantially true contour 37 for the profile 35 is 
defined. The sharpness of the true contour 3 7 is in 
contrast to the smeared contour that without de- 
convolution would have been obtained by virtue of the 
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spread (represented by the arrows 2 6) of the relevant 
image ERFs 25. 

A small portion of the contour 37 is shown enlarged in 
5 Figure 4 and is defined as a best-fit line between 

optimum edge-positions derived respectively from x- and 
y-profiles; the x-profile positions are indicated by 
black dots and the y-profile positions by white dots. 
The closeness to one another of corresponding positions 

10 in the x- and y-profiles is indicative of the accuracy to 
sub-pixel level achieved. The smear-out that would have 
been manifest in the image-profile contour if the de- 
convolution technique were not used, would have extended 
throughout the space bounded by the dashed lines ^3 8 and 

15 39; these boundaries are indicated in dashed line in the 
representation of profile 3 5 in Figure 2. 

The accurate definition of the image contour 37 derived 
in the display 8 allows correspondingly accurate 

20 determination in calculation 9 of the area within that 
contour; the volume involved can also be derived from 
successive slices scanned. The determination of area and 
volume is especially useful for diagnostic and accurate 
assessment 10 of the size of a tumour or lesion before 

25 and after treatment. It is similarly useful for 
assessment of arterial dimensions in angiography. 

Moreover, the accurate definition of the image contour 3 7 
derived in the display 8, is particularly useful for 

30 segmenting anatomical structures for diagnostic and 
treatment planning 11. Furthermore, the ratio of 
intensities of two scans are derived by processing 12 to 
derive values of relaxation times T, and T 2/ and of proton 
density. The values are then represented in display 13 

35 within the boundary of the image contour, utilising 

standardisation data derived from a couch-top 14 used 
within the scanner 1. The couch-top JL4, which also 
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provides landmarks for determining position coordinates, 
has the form shown in Figures 5 and 6, and will now be 
described. 

Referring to Figures 5 and 6, the couch-top 40, which is 
of polystyrene bubble foam, has the form of a flat slab, 
and is supported on a curved foam base 41 that fits to 
the scanner table of the MR installation. Two zig-zag 
runs of tubing 4 2 are embedded within the top 4 0 to 
extend along the length of the couch either side of 
centrally-embedded straight tubing 43. 

The tubing 42 of each zig-zag run is of double-bore 
rectangular cross-section, whereas the tubing 4 3 is of 
single-bore cross-section. The five bores defined by the 
array of tubing 42 and 43 may be filled respectively with 
the five MR solutions S 0 to S 4 of Table I, for 
standardisation and calibration purposes. The four 
MnCl 2 .4H 2 0 solutions, S, to S 4 , cover the full range of 
values of Tj and T 2 for anatomical tissues, and the fifth 
solution, S 0 , of CuS0 4 .5H 2 0, is nominally equivalent to 
"loosely bound water". 



Table I 



Solution 




at 0.5T 


T 2 

at 0*5T 




1.25 g/1 CuS0 4 .5H 2 O 


200 ms 


2 00 ms 


s, 


3.41X10 16 Mn +2 ions/ml 


84 0 ms 


300 ms 


S 2 


1.15xl0 17 Mn +2 ions/ml 


44 0 ms 


120 ms 


s 3 


2.30X10 17 Mn +2 ions/ml 


2 50 ms 


6 0 ms 


s 4 


4.37X10 17 Mn +2 ions/ml 


150 ms 


3 0 ms 
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The tissue types revealed by the T,, T 2 and proton density 
values in display 13 are determined by processing 15 from 
look-up tables, and tissue densities are assigned within 
the image-contour boundaries in display 16. The images 
of display 16 are furthermore corrected empirically for 
geometry distortion in accordance with data supplied from 
memory 17. The data stored in memory 17 is derived using 
a set of drum phantoms of the spider-web form, and 
correction for geometry distortion is realistically 
effective principally because of the accuracy and spatial 
resolution with which image-contours are defined. 

The tissue types assigned to the corrected images are 
utilised in display 16 through the Bulk Heterogeneity 
Correction method described by Richard A , Geise et al , 
Radiology, 124:133-141, July, 1977, to establish for each 
image a normalised tissue density value; the up-datable 
look-up table for this is stored in memory 18. 
Accordingly, the display 16 when used at step 19 in 
conjunction with the positional datums derived from the 
couch 14, has all the tissue contours accurately mapped 
out with their respective tissue densities and locations. 
This establishes an accurate and readily-usable, stand- 
alone basis for diagnosis and treatment planning, and 
enables a true three-dimensional assessment and plan to 
be made when both orthogonal and oblique MR images are 
involved. 

Although the method and system of the invention have been 
described above in the medical context they are 
applicable more widely than this, for example, in 
engineering, in physical science and in the field of 
instrumentation generally. Moreover, the method and 
system is not limited to MR imaging, but may be utilised 
where other forms of imaging are involved. The steps and 
structure represented in Figure 1 by x boxes' l to 11 are 
just as applicable to computer assisted tomography (CT) , 
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as they are to MR imaging. Other forms of imaging to 
which the invention is applicable include X-ray 
radiography, film- or print-image transformation to 
digital form, digital X-ray f luorography , ultra-sound 
5 imaging, nuclear medicine, positron emission tomography 
(PET) and other camera or imaging. The technique is 
particularly suitable for use in X-ray digital 
f luorography, in which small structures under study are 
highlighted by injection of contrast liquids; the small 
10 structures may also be isolated from surrounding 
interfering effects by using an image subtraction 
technique. 

The inherent resolutions of X-ray radiography, 
15 ultra-sound imaging, nuclear medicine, and PET scanning 
are relatively low, and some are used for real-time 
study. Only the individual still frame or hard-copy 
images may be re-processed. 

2 0 In the context of engineering, physical science and the 
field of instrumentation, the invention is applicable to 
one-dimensional imaging as used, for example, in regard 
to bar-code patterns, the spectrum of DNA analysis, iris 
patterns of eyes (for example, for identification 

25 purposes in commercial banking) , finger-print 

identification, and emission spectroscopy. The invention 
is also applicable to two-dimensional imaging, for 
example, in relation to images obtained by satellite or 
pattern recognition, or from a surveillance camera or 

30 during laboratory experimentation. As a general matter, 
the invention is applicable where there needs to be 
accurate determination of the edge position in an image 
versus the true object-edge position, for the purpose, 
for example, of measurement of the positional 

35 displacement between object and image, distortion 
correction and manufacturing control. 
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As a further example of application of the present 
invention, a method and system that uses CT and MR 
imaging in conjunction with one another, will now be 
described . 

The major contribution to the magnetic-resonance (MR) 
signal comes from the abundant protons content of water 
molecules and protein. It is a quantum process at the 
Larmor frequency according to the magnetic field in use. 
The x T,-weighted' and x T 2 -weighted' MR signals from 
protons provide contrast numbers that are relative in 
scale, whereas in CT, the X-ray absorption is a 
polychromatic attenuation process affected by the 
electron densities of all the atoms presented within the 
X-ray beam. There is no equation to correlate the CT 
number (or the linear attenuation coefficient, electron 
density, or tissue density) with the MR-contrast numbers; 
no direct calibration between the two types of signal is 
possible. This lack of correlation is confirmed by 
consideration of bone and air which are at opposite ends 
of the CT contrast (absolute) scale using water as the 
base-line reference, but which are at the same end of the 
MR-image contrast (relative) scale owing to their common 
low proton-population. 

The lack of correlation between the CT and MR signals 
acts against their use in combination for imaging 
purposets, but the present invention provides a method and 
system by which the advantages of each may be utilised to 
improve image resolution and contrast information. 

In the latter regard, CT provides a high spatial 
resolution but only in regard to view normal to the 
transverse slice. Resolution for all reconstructed 
non-transverse planes is poor owing to the need to use 
elongate voxels to improve signal-to-noise ratio. Also, 
the partial-volume effect of using elongate voxels may 
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give rise to detection errors at the thin edge of a 
contrast profile of a lesion. MR, on the other hand, can 
give the same high degree of spatial resolution viewed in 
the normal direction to any image-slice plane, and can 
also provide isotropic resolution with cubic-voxel volume 
imaging. 

To this end, multiple-slice transverse CT scans are 
collected across a section of the volume of interest in a 
patient or other subject. Corresponding multiple-slice 
transverse MR scans of the same volume are also 
collected. The slice thickness of the latter scans may 
be one half, or smaller, of the thickness of the CT 
slices, and may be collected two-dimensionally or three- 
dimensionally . The patient or other object scanned is ^ 
constrained throughout on a couch that provides landmarks 
with respect to a coordinate reference arrangement on the 
couch-top. This is to ensure the reproducibility of, for 
example, anatomical positions and features to the first 
order accuracy for the corresponding CT and MR scans, and 
possibly for radiation treatment to be made. 

The respective transverse planes of the CT and MR images 
are processed individually in the method and are matched 
with,, one another in a de-convoluted space for the CT and 
MR images. The two sets of de-convoluted maps are then 
merged together to a second order of accuracy in order 
that the CT numbers may be transferred over to replace 
the corresponding MR contrast numbers. Once this has 
been achieved, non-transverse (or oblique) planes can be 
obtained from the two-dimensional MR images, or from the 
re-arrangement of the corresponding voxels of the three- 
dimensional volume images; where two-dimensional MR is 
used, a further step of contrast transformation may be 
required. 
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The respective transverse planes of CT and MR images are 
processed individually by using x boundary' or 
x finger-print' matching techniques in a de-convolution 
space for the CT & MR images. In these transverse image 
planes, in particular when used medically, the 
skin-contour features along the sides of the patient may 
be best used for second-order alignment and matching 
purposes, as they are less affected by patient-movements. 
Transitional error may be readily corrected with respect 
to the coordinate positions of a rectangular tubing 
system embedded in the couch-top (for example that 
described above with reference to Figures 5 and 6) . The 
processed data may then be used for a N diagnostic and 
statistics software package' of CT image versus MR image 
for their exactly corresponding transverse slice (s), and 
an associated % statistical package' for accurate 
computation of the 'true' area, and then the 'true' 
volume, of a lesion or tissue profile or contour. 

The two sets of de-convoluted maps may also be merged 
together to a second order of accuracy in order that the 
CT numbers may be transferred over to replace the 
corresponding MR contrast numbers. Once this has been 
done, non-transverse (or oblique) planes are obtained 
from the two-dimensional image or from the re-arrangement 
of the corresponding voxels of the three-dimensional 
volume; in the two-dimensional MR case, a further step of 
contrast transformation is required. The transferred 
contrast data may then be used in a three-dimensional 
radiotherapy treatment planning software package for an 
in-plane, oblique-image pseudo-three-dimensional approach 
using MR images. 

Software required for de-convolution processing of image 
data according to the invention may be implemented in 
conjunction with a least-squares curve-fitting method. A 
mathematical model of the method, fro*n which the required 



WO 00/22573 



PCT/GB99/034I7 



1.5 

software can be readily developed, will now be given in 
relation to Figure 7 which indicates the convolution of 
an LSF L(x) with a step function SF having a step-down 
edge SDE , to produce an image ERF E(x) . Considering the 
values of E(x) resulting from the convolution at points 
x=l to x=n of SF: 



at x=l : 



£i(x) Ax = E(x) x 



at x=2 



^L(x)Ax - ]Tl(x)Ax = E(x) 



a-1 



at x=3 : 



J2L(x) Ax - £ L(x) Ax = E(x) 3 



a-2 



until at x=n: 



J2L(x)Ax -^L(x)Ax = E(x) n 



For de-convolution, that is to say, for the ERF E(x) to 
manif est the LSF L(x) : 



E(x) 1 - E(x) 2 = JTL(x)Ax = L(x) (a . u . a Ax 

a-l 

{E(X) X - F(X) 2 )/AX = L(X) (a -i)-a 

thus: 

{E(X) 2 - E(x)J/Ax = L(X) ( a- 2 )-<a-l> 

leading to 

iE(x) (n . 1} - F(x) n )/Ax = L(x) {nml) _ n 
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Thus, the same shape of LSF is recovered by the 
de-convolution process independently of the sense of the 
ERF; it does not matter whether a roll-up or a roll-down 
ERF of a positive or a negative contrast profile, is 
involved. This is an important property as a positive 
contrast contour will have roll-up (from low to high) 
ERFs at both ends whereas a negative contrast contour 
will have roll-down (from high to low) ERFs at both ends. 

In the practical implementation , the LSF is derived from 
the ERF by de-convolution using a running filter. This 
enhances the accuracy of the method in overcoming the 
problem of noise that affects the digitised pixel values 
of the image. Use is made of a least-squares fitting of 
a set of data points that advances along the whole length 
of the function from one point to another. 

Assuming that y = a 0 + a,x + a 2 x 2 represents the ERF curve 
a five-point or seven-point fit is used, and the normal 
equation becomes: 



Ei 


Ex s 


Ex! 




aO 




Ex, 


Ex! 


Ex! 




al = 


= Ex iYi 


Ex! 


Exj 


Exj 




a2 





where: i is 1+n, ... 5+n; ... until m-5, ... m 

or 

i is l+n, . . . 7+n; . . . until iii-7, ... m 
n is 0, 1, 2 , ... ; and 

1, ... m is the span of the ERF profile. 
The solution may be derived from either: 



aO 


El 


Ex 


Ex 2 


-1 


Ey- 


al ■- 


- Ex 


Ex 2 


Ex 3 




Exy 


a2 


Ex 2 


Ex 3 


Ex 4 




Ex 2 y 
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or: 



aO 



Ey Ex Ex 2 
Exy Ex 2 Ex 3 
Ex 2 y Ex 3 Ex 4 



al 

El Ey Ex 2 
Ex Exy Ex 3 
Ex 2 Ex 2 y Ex 4 



a2 



Ey 

Exy 
Ex 2 Ex 3 Ex 2 y 



Ei Ex 

Ex Ex 2 



Ei Ex Ex : 
Ex Ex 2 Ex 
Ex 2 Ex 3 Ex 4 



For both of these equations to be valid :- 



El Ex Ex 2 
Ex Ex 2 Ex 3 
Ex 2 Ex 3 Ex 4 



* 0 



5 The gradient at dy/dx (3+n) or dy/dx (4+n) can then be derived 
and plotted against x for the LSF profile. 

The graph of dy/dx against x gives the LSF profile. The 
peak of this profile is located centrally of the mid- 
10 points of the ascending and descending limbs of the 

graph. These points define the extremes of the full- 
width-half-maximum (FWHM) of the profile and the mid- 
point of this is determined with an accuracy of sub-pixel 
level owing to the * average' effect. 

15 

The point spread function (PSF) is the two dimensional 
profile which may be derived, in practice, from the two 
corresponding LSFs orthogonal to one another within an 
image plane. The peak position of the PSF profile is, 
20 therefore, from the x mean' or ^cross-over' of the two 
peaks or the two LSF profiles. The PSF is obtained in 
practice from two orthogonal axes in a two-dimensional 
plane . 

25 The generation of the LSF (or PSF) is, after phase- 
reversal correction, independent of the roll-up or 
roll-down nature of ERFs at the edges of the contrast 
contour* In other words, it is independent of the sense 
and the absolute value of the contrast* numbers within the 
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contour. The peak position of the LSF (or PSF) is the 
central half-way (50%) point of the roll-up or roll-down 
ERF which is the optimum position for true-edge 
definition. 



